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q ' Abstract. Gonzalez etal. (2003) have reported the discovery of an anomalous radiation component 

, from f» 1 - 200 MeV in GRB 941017. This component varies independently of and contains £ 3 x 

■ the energy found in the prompt ~ 50 keV - 1 MeV radiation component that is well described by the 
relativistic synchrotron-shock model. Acceleration of hadrons to very high energies by GRBs could 

ON , give rise to a separate emission component. Two models, both involving acceleration of ultra-high 

energy cosmic rays with subsequent photomeson interactions, are considered. The first involves 
a pair-photon cascade initiated by photohadronic processes in the GRB blast wave. Calculations 
indicate that the cascade produces a spectrum that is too soft to explain the observations. A second 
ly-^ ' model is proposed where photopion interactions in the GRB blast-wave shell give rise to an 

^vq escaping collimated neutron beam. The outflowing neutrons undergo further photopion interactions 

■ to produce a beam of hyper-relativistic electrons that can lose most of their energy during a fraction 
7-H ' of a gyroperiod in the Gauss-strength magnetic fields found in the circumburst medium. This 

, secondary electron beam produces a hard synchrotron radiation spectrum that could explain the 

■ anomalous component in GRB 941017. 

2a: 

Of 

6 INTRODUCTION 

a: 

Based on joint analysis of BATSE LAD and EGRET TASC data, Gonzalez et al. (2003) 
recently reported the detection of an anomalous MeV emission component in the spec- 
trum of GRB 941017 that decays more slowly than the prompt emission detected with 
the BATSE LAD in the « 50 - 300 keV range. The multi-MeV component lasts for 
^ 200 seconds (the % duration of the lower-energy component is 77 sec), and is de- 
tected with the BASTE LAD near 1 MeV and with the EGRET TASC between « 1 and 
200 MeV The spectrum is very hard, with a photon number flux 0(e) oc £~\ where 
£ = hv/m e c 2 is the observed dimensionless photon energy. 

This component is not predicted or easily explained within the standard leptonic 
model for GRB blast waves, though it possibly could be related to Comptonization 
of reverse-shock emission by the forward shock electrons 01, including self-absorbed 
reverse-shock optical synchrotron radiation [11]. Another possibility is that hadronic 
acceleration in GRB blast waves could produce this component. 

We consider two models involving acceleration of hadrons at the relativistic shocks of 
GRBs. In the first model, the anomalous radiation component would be explained by the 
pair-photon cascade radiation initiated by photohadronic processes between high-energy 
hadrons accelerated in the GRB blast wave and the internal synchrotron radiation field. 
We find, however, that the cascade radiation spectrum is too soft to be able to explain 
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FIGURE 1. Photon energy fluence from an electromagnetic cascade initiated by photopion secondaries 
in a model GRB, with parameters given in the text. Five generations of Compton (heavy curves) and 
synchrotron (light curves) are shown. The first through fifth generations are given by solid, dashed, dot- 
dashed, dot-triple-dashed, and dotted curves, respectively. The total cascade radiation spectrum is given 
by the upper bold dotted curve. 



the anomalous emission component of GRB 941017. In the second model, photomeson 
interactions in the relativistic blast wave produce a beam of ultra-high energy neutrons, 
as proposed for blazar jets |2|]. Photopion production of these neutrons with photons 
outside the blast wave produce a directed hyper-relativistic electron-positron beam in 
the process of charged pion decay and the conversion of high-energy photons from n° 
decay. These energetic leptons produce a synchrotron spectrum in the radiation reaction- 
limited regime extending to ^ GeV energies, with properties in the 1 - 200 MeV range 
similar to that measured from GRB 941017. If this model is correct, detection of this 
component therefore gives important indirect evidence for the acceleration of ultrahigh 
energy cosmic rays in GRB blast waves. 



HADRONIC MODELS FOR GRB 941017 

We assume that efficient proton acceleration to ultra-high energies takes place in GRB 
blast waves, as required in models where GRBs accelerate high-energy cosmic rays 
Photopion interactions of these protons with internal synchrotron photons 
and with photons from an external radiation field create neutral particles and charged 
pions that decay and initiate an electromagnetic cascade within the GRB blastwave. 

One possible model is to attribute the anomalous component to radiation of the 
cascade produced within the blast wave. Fig. 1 shows the hadron-initiated cascade 
radiation for a model GRB at redshift z = 1, with hard X-ray fluence 4> /of = 3 x 
10~ 4 erg cm~ 2 , a light curve of 100 second duration divided into 50 pulses of 1 second 
each, and Doppler factor 8 = 100 (see ||4i |2J] for more details about the model). The 
total amount of accelerated proton energy E' = 4ndj<i> tot 8~ 3 (l +z) _1 is injected into 
the comoving frame of the GRB blast wave. 



The various generations of synchrotron and Compton radiation initiated by the cas- 
cade are shown in Fig. 1, along with the total radiation spectrum. As can be seen, the 
isotropic cascade radiation approaches the spectrum of an electron distribution cooling 
by synchrotron losses, that is, a spectral index between —1.5 and —2. This is too soft 
to explain the observations of GRB 941017 with a spectral index = — 1. Moreover, the 
isotropic cascade radiation should decay at the same rate as the synchrotron radiation. 
Thus the radiation from an electromagnetic cascade in the GRB blast wave cannot ex- 
plain the anomalous component in GRB 941017. 

We propose a second possibility based on our neutral beam model J2]. Ultra- 
relativistic neutrons formed in the reaction p + y — > n + n° are not confined by the 
magnetic field in the GRB blastwave shell and flow out to create an energetic neutron 
beam. These neutrons are subject to further photopion processes with photons in the 
surrounding medium to form charged and neutral pions. In the Gauss-type magnetic 
fields surrounding GRB sources that we assume here, charged n and \x at energies 
;$ 10 18 eV decay rather than lose energy through synchrotron emission JH. The 
charged pions decay into ultrarelativistic electrons and neutrinos, whereas the decay 
of n° produces two y rays that are promptly converted into electron-positron pairs on 
this same external radiation field. These energetic electrons (including positrons) are 
initially produced in the direction of the GRB jet. 

The spectra of secondary electrons created by the neutron beam displays a sharp 
cutoff at energies < 10 14 eV as a consequence of the high threshold for photomeson 
interactions (see Ref. flj] for more details). Electrons with Lorentz factor y lose en- 
ergy through synchrotron radiation in an ordered magnetic field with strength B at the 
rate —dy/dt = OjB 2 y 2 sin 2 y// (4zm e c), where \]/ is the electron pitch angle. The corre- 
sponding synchrotron energy-loss time scale t syn = y/\dy/dt\. The gyration frequency 
(Ob = eB/m e cy, and is independent of pitch angle. When (OBt syn <C 1, the electron loses 
almost all of its energy into synchrotron radiation in a time less than the gyroperiod. We 
use the term "hyper-relativistic" to refer to electrons in this radiation-reaction regime of 
synchrotron emission ful fill. The mean energy of synchrotron photons from electrons 
that enter the hyper-relativistic regime is independent of magnetic field (e.g., J9J). 

Electrons which cool before rotating by an angle equal to the jet opening angle 0j will 
produce synchrotron photons that are primarily directed within the jet opening angle. 
This condition is defined by (Qstsyn < Qj, which applies to electrons with y ^ Yhr(Qj) = 



4ne/(djG T Bsm 2 iff) = 3 x lO 7 /[sini// v /(0 / /O.l)5(G)]. Here we have taken a typical 
jet opening half-angle Qj = 0. 1 because this is the average value implied by the analysis 
leading to the standard energy reservoir result of GRBs [5]. Lower-energy electrons with 
7 < Yhr{®j) radiate their energy over a much larger solid angle and longer time. 

The characteristic synchrotron photon energy Sy = m e c 2 £ radiated by electrons which 
lose their energy within the jet opening angle Qj is 

heBsmy 500 
] m e c (1+z) (e//0.1)[(l+z)/2]sinvf ' K) 



Hyper-relativistic electrons with y > % r (Qj) rapidly lose energy through synchrotron 
losses and deposit all of their energy along the direction of the jet. Electrons at lower 



energies are deflected to angles 6 > Oj, and their emission is not seen by an on-axis 
observer. Hence the distribution of electrons along the jet direction always has an 
effective low-energy cutoff at % r (0y). The production spectrum of the electrons can 
have an intrinsic cutoff y co due either to the low-energy cutoff in the escaping neutron 
spectrum, or to the neutron-induced photomeson secondary spectrum outside the GRB 
blast wave. If p = y C o/Yhr(9j) > 1» then the observed synchrotron spectrum is a power 
law with —1.5 index for Sj ^ Sy £ p Sj, and a photon spectrum with the same spectral 
index as the accelerated protons and escaping neutrons at photon energies Sy ^ p 2 Sj 
ifioll . If p < 1, then the observed photon spectrum at Sy £ Sj has the same spectral index 
as the primary neutrons. 

At photon energies Sy -C Sj, the observed spectrum is produced by the same hyper- 
relativistic electrons with y^ yh r (0j), but at energies £ well below the peak energy 
3y 2 £ B , where £ B = B/B cr and B cr = 4.41 x 10 13 G is the critical magnetic field. We 
now derive this spectrum. 

The differential energy radiated per dimensionless energy interval de per differential 
solid angle element dQ. in the direction 6 with respect to the direction of an electron 
moving with Lorentz factor y is given by 

dE e ~ ■(^) 2 (i + r 2 2 ) 2 (Aii+A±), (2) 



dedQ. 3 n 2 Ac Y £ b' 

where Ac = k/m e c = 3.86 x 10~ n cm is the electron Compton wavelength, and An = 
Ky 3 (<j;) and A_l = (yO) 2 K 2 ^ 3 (£,)/[l + (yd) 2 } are factors for radiation polarized parallel 
and perpendicular to the projection of the magnetic field direction on the plane of the 
sky defined by the observer's direction [8]. The factor £ = e/e, where £ = 3£b7 2 /(1 + 
y 2 9 2 ) 3 / 2 , andK„(x) is a modified Bessel function of the second kind, with asymptotes 
K n (x) — > jT(n)(2/x) n in the limit x<l, and K n (x) — > ^7t/2x exp(— x) in the limit 
x 3> 1. The condition £ <C 1 corresponds to £ £ where K n (£) are in their power- 
law asymptotes, and £ ^ 1 or £ £ £ is where K„(£) are in exponential decline. The 
characteristic energy £ approaches 3£gy 2 when yO <^1, and £ declines with according 
to the relation £ = 3e B y 2 /(yO) 3 when yd > 1. When £ < £, then An > A_l and 

dE/dedQ = (dE syn /dedQ.) ~ 3 l / 3 (l.07e/n) 2 (ye/e B ) 2 / 3 /Ac £ 2/3 - For a fixed value 
of £, this emissivity exponentially cuts off when £ > 3e B /y0 3 , or when £ max = 
(3£ B /ye) 1 / 3 . 

The synchrotron emission spectrum in the limit Sy <C Sj, integrated over solid angle, 
is simply given 

j£ Jo K dedQ J nlc 

This differs from the energy index +1/3 for the electrons in the classical regime aver- 
aged over a complete orbit because, in this case, dO s'mO — > dO sini// in the integration 
ineq. © id. 



NEUTRON BEAMS FROM GRBS 



In our model, ultrarelativistic protons undergo photomeson interactions with the internal 
and external radiation photons [4], producing a beam of outflowing neutrons. Subsequent 
interactions of these neutrons with the external radiation field generates a beam of hyper- 
relativistic electrons. (The external photons could be due to pre- supernova thermal or 
GRB synchrotron radiations scattered by the stellar wind in the collapsar model, or by 
plerionic emission in the SA model.) The synchrotron radiation has a low-energy cutoff 
at Sy < Sj ~ 500 MeV - 1 GeV, with a specific characteristic power-law number spectral 
index equal to —1, as observed for GRB 941017 [6]. 

The ps 200 s decay time of the anomalous emission can be explained by the emis- 
sion of hyper-relativistic electrons from the edges of a jet blastwave at distances R ~ 
1O 15 (0//O.1) _2 [(1 +z)/2] _1 cm, implying significant opacity to photomeson processes 
due to external photons at these distances. The differing GRB external radiation and 
density environments which determine the intensity of target photons could account for 
the unusual spectrum of GRB 941017. The neutron beam can carry up to ~ 50% of the 
entire energy of injected protons J2J], which can be further reprocessed about equally into 
neutrinos and hyper-relativistic electrons. This nevertheless requires a large nonthermal 
hadron-to-lepton load in the GRB blast wave in order for the fluence of the anomalous 
component to be comparable with that of the normal component. This large hadron load 
is also consistent with predictions linking cosmic rays to GRB sources IU6I1 . 

Our model predicts that most of the energy of hyper-relativistic electrons will be 
observed at ^ 500 MeV - GeV energies as a delayed radiation component, well-suited 
for observations with GLAST. GRBs with anomalous y-ray emission components should 
also be bright neutrino sources detectable with IceCube. 
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